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Abstract 



a.: 

■ We present Version 3.0 of the GLoBES ("General Long Baseline Experiment Sim- 

ulator") software, which is a simulation tool for short- and long-baseline neutrino 
oscillation experiments. As a new feature, GLoBES 3.0 allows for user-defined system- 
atical errors, which can also be used to simulate experiments with multiple discrete 



sources and detectors. In addition, the combination with external information, such 
as from different experiment classes, is simplified. As far as the probability calcula- 
tion is concerned, GLoBES now provides an interface for the inclusion of non-standard 
physics without re-compilation of the software. The set of experiment prototypes com- 
ing with GLoBES has been updated. For example, built-in fluxes are now provided for 
the simulation of beta beams. 
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1 Introduction 



Neutrino oscillations are now established as the leading flavor transition mechanism for 
neutrinos in a long history of many experiments, see e.g. Ref. [1] and references therein. 
Future facilities, using accelerator-based neutrino beams or nuclear reactors as neutrino 
sources, are proposed for precision measurements of the neutrino oscillation parameters. In 
the simulation of these experiments, the presence of multiple solutions which are intrinsic to 
neutrino oscillation probabilities [2-5] affect the performance. Thus, optimization strategies 
are required which maximally exploit complementarity between experiments. The GLoBES 
software package [6] is a modern experiment simulation and analysis tool for a highly ac- 
curate beam and detector simulation. In addition, it provides powerful means to analyze 
correlations and degeneracies, especially for the combination of several experiments. Com- 
pared to a Monte Carlo simulation, which yields a different result in each run, it simulates 
the performance of the average experiment (see Ref. [7] for a discussion of the meaning of 
"average"). The advantage of such an average prediction is a tremendous performance gain, 
which can be used for systematical parameter space scans. In addition, it simplifies the 
direct comparison of experiments. 

The GLoBES software has, in the past, been used for many studies, some of them are re- 
ferred to in the paper. However, recent developments have indicated that extensions and 
improvements are necessary. In this work, we present the most important new features and 
changes in GLoBES 3.0. For experimentalists, GLoBES now allows for the implementation 
of arbitrary systematics, which can also be used for the simulation of multi-detector exper- 
iments. In addition, several extensions have been included in AEDL ("Abstract Experiment 
Description Language"), such as built-in beta beam fluxes and the possibility to use lists as 
variables. For phenomenologists, user-defined priors provide a flexible interface to include 
external information, such as from different experiments. In addition, GLoBES now can be 
used for the simulation of non-standard physics without re-compiling the GLoBES software. 

2 Concept of GLoBES 

GLoBES ("General Long Baseline Experiment Simulator") is a flexible software tool to sim- 
ulate and analyze neutrino oscillation short- and long-baseline experiments using a complete 
three-flavor description. On the one hand, it contains a comprehensive abstract experiment 
definition language (AEDL), which allows to describe most classes of long baseline and reac- 
tor experiments at an abstract level. On the other hand, it provides a C-library to process 
the experiment information in order to obtain oscillation probabilities, rate vectors, and 
Ax 2 -values (c/., Fig. CD). In addition, it provides a binary program to test experiment def- 
initions very quickly, before they are used by the application software. Currently, GLoBES 
is available for GNU/Linux. 

GLoBES allows to simulate experiments with stationary neutrino point sources, where each 
experiment is assumed to have only one neutrino source. Such experiments are neutrino 
beam experiments and reactor experiments. Geometrical effects of a continuous source 
distribution, such as in the sun or the atmosphere, can not be described. In addition, 
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Figure 1: General concept of the GLoBES package. 



sources with a physically significant time dependence can not be studied, such as supernovae. 
However, in GLoBES 3.0 and higher, new flexibility is introduced by the concept of user- 
defined systematics. This new feature allows the cross-definition of systematical errors over 
different experiments. In principle, this mechanism can be used for the combination of 
several discrete sources and one detector, or several detectors and one source, or several 
sources and several detectors. In addition, already implemented concepts of GLoBES have 
been used for indirect simulations of geometrical effects. For example, the mapping of the 
detector location on the neutrino energy has been simulated in Ref. [8] by the use of variable 
bin widths. 

On the experiment definition side, either built-in neutrino fluxes (e.g., neutrino factory, 
beta beam) or arbitrary, user-defined fluxes can be used. Similarly, arbitrary cross sections, 
energy dependent efficiencies, energy resolution functions as well as the considered oscillation 
channels, backgrounds, and many other properties can be specified. For the systematics, 
energy normalization and calibration errors can be simulated, or the systematics can be 
completely user-defined. Note that energy ranges and windows and bin widths can be 
(almost) arbitrarily chosen, including bins of different widths. Together with GLoBES comes 
a number of pre-defined experiments in order to demonstrate the capabilities of GLoBES 
and to provide prototypes for new experiments. In addition, they can be used to test new 
physics ideas with complete experiment simulations. 

With the C-library, one can extract the Ax 2 for all defined oscillation channels for an 
experiment or any combination of experiments. Of course, also low-level information, such 
as oscillation probabilities or event rates, can be obtained. GLoBES includes the simulation 
of neutrino oscillations in matter with arbitrary matter density profiles. In addition, it allows 
to simulate the matter density uncertainty (see, e.g., Refs. [9,10]) and to extract the precision 
on the matter density (see, e.g., Refs. [11,12]). As one of the most advanced features of 
GLoBES, it provides the technology to project the A% 2 , which is a function of all oscillation 
parameters, onto any subspace of parameters by local minimization. This approach allows 



2 



the inclusion of multi-parameter-correlations, where external constraints (e.g., on the solar 
parameters) can be imposed, too. Applications of the projection mechanism include the 
projections onto the sin 2 2^13-axis and the sin 2 2#i3-(!>cp-plane. In addition, all oscillation 
parameters can be kept free to numerically localize degenerate solutions. 

3 Oscillation probabilities and the simulation of non-standard 
physics 

The probability calculation in GLoBES is based on the diagonalization of the Hamiltonian 
in layers of constant matter density using the standard three flavor scenario of neutrino 
oscillations. In the flavor base, we have 

, ( \ / ±V2G F n e \ 

H(n e ) = —U Am 2 , \rf + (1) 

Ztj \ Am 2 x / V / 

with U being the mixing matrix (described by four parameters #12, #13, #23, and 5cp), and 
n e the constant electron density in the respective matter density layer. The electron density 
is related to the matter density p by n e ~ 0.5 p/m^ with being the nucleon mass. The 
sign in the second term depends on whether one uses neutrinos or antineutrinos. For k 
matter density layers with densities n k and thicknesses Xj, the oscillation probability then 
evaluates to 

P(v a -> up) = \(v p \£(n k ,x k ) . . .£(n\,xi)\v a )\ 2 (2) 
with the evolution operators (see, e.g., Ref. [13]) 

e -im1x k 
£(n k e ,x k ) = exp(-iH(n k )x k ) = — U(n k e ) ( e~ im ^ | U(n^ . (3) 



e -im^x k 



Here the Hamiltonian 7i(n k ) is diagonalized by the unitary mixing matrix in matter U(n k ) 
with the eigenvalues m k , m k , and m k . The probabilities in Eq. (j2J) are used for the event 
rate computation, which is described in greater detail in Ref. [6]. In the next section, we 
will demonstrate how systematics is implemented. 

GLoBES 3.0 allows the modification of this Hamiltonian and the whole probability engine 
without re-compilation of the GLoBES software. This is implemented by using pointers to 
the functions of the standard probability engine, which can be changed by registering a 
different probability engine. For example, a user may just copy the standard probability 
engine of GLoBES, modify it, and register it. An important prerequisite is the ability 
to handle more than the standard six oscillation parameters (plus matter density) A = 
(6*12 , 6*i3 , 6*23 , <5cp , Ar/4 , Am 2 x ,p), i.e., 



> '12 



6*13 , #23 , <$cp , Am 2 x , Am 2 31 , p, r) X , . . . , r} n ) (4) 



with n non-standard parameters m. These parameters can be accessed in the usual way, 
see Ref. [15]. For example, non-standard Hamiltonian effects in the e-r-sector may be 
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Figure 2: Fit in the sin 2 2#i3-<7E-plane for the simulated values sin 2 26*i3 = erg = (figure similar to 
Ref. [14]). Here oe is an additional non-standard parameter describing the wave packet width of the neutrino 
oscillation wave packet. If this parameter is kept free (vertical axis), the sin 2 2 9 13 sensitivity (horizontal 
axis) will be affected. Note that the standard sin 2 26*13 sensitivity is given by the curves for erg = 0. The 
figure is computed with example6.c coming with the GLoBES 3.0 distribution (and higher). 

introduced by modifying Eq. ([TJ) to 

/ e eT \ 

H = H{n e ) + (5) 

\4 / 

In this case, there are two more real parameters, such as the absolute value and phase of 
e er . In the literature, the non-standard physics feature in GLoBES as experimental feature 
has been used for damping effects (such as neutrino decoherence, decay, etc.) in Ref. [14], 
Hamiltonian-level effects (such as non-standard matter effects) in Ref. [16], and mass- varying 
neutrinos in Ref. [17] by implementing an environment dependence of neutrino mass. We 
show in Fig. [2] a result from Ref. [14] provided as example6.c in the GLoBES distribution. 
In this case, the non-standard physics is a loss of coherence in neutrino oscillations described 
by an intrinsic width of the mass eigenstate wave packets oe (in energy space). Compared 
to Eq. (j3J), the underlying physics takes place at the probability level, i.e., Eq. ([2]), and 
it can be easily implemented analytically (for details, see Ref. [14]). In example6.c, an 
analytical implementation of this physics is used for simplicity. 

4 Systematics implementation and user-defined systematics 

GLoBES supports four types of systematical errors by default: Signal and background nor- 
malization errors, as well as signal and background tilt or energy calibration errors. The tilt 
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(T) is implemented as a linear distortion of the spectrum around the center, and the energy 
calibration (C) as a distortion (stretching) of the reconstructed energy scale. For the im- 
plementation of these systematical errors, the "pull method" is used [18], which introduces 
so-called nuisance parameters Q. For example, for the signal and background normalization 
errors, the respective event rates s, and 6j in each bin i are multiplied with0 

Si(Ci) == (1 + Ci) • Si , b t (( 2 ) := (1 + Ca) ■ h , (6) 

i.e., the rates are scaled by these parameters. The systematics x 2 is then minimized over 
these parameters: 

Xpuii(A) := mm ( X 2 (A, Ci, . • • , Q + E 4 ) ■ (7) 

{U V 3=1 °v / 

Here x 2 (A, Ci; • • • > Cfc) i s the usual Possonian x 2 depending on the neutrino oscillation param- 
eters A and the nuisance parameters Q (which, for instance, scale the signal or background 
rates). In addition, Gaussian penalties C]/^ are added, where <tq corresponds to the ac- 
tual systematical error. The core part of Eq. (JTj) is the definition of \ 2 . For example, for a 
background-free measurement and a signal normalization error only, it reads 

# of bins , n \ 2 

X 2 (A,a)= 2(^(l + a)T t -O t + Oaog (1 + ; )T J (8) 

where Oi are the observed rates (corresponding to the data or true values), and T, are the 
theoretical (fit) rates. Let us assume that we have two detectors, such as for a reactor 
experiment with a near and far detector. Using standard systematics in GLoBES, one has 

xJuii(A) = mm ^ x 2 (A,a) + ^ + nun (xt(\b) + ^ . (9) 

In this case, the normalization between the two detectors will be uncorrelated. 

Using GLoBES 3.0 and higher, one can modify the standard x 2 function in GLoBES which 
corresponds to Eq. (jSJ). In addition, the systematics concept of "error dimensions" is re- 
placed by a concept directly related to this % 2 function. As an example, consider a reactor 
experiment with near (N) and far (F) detectors which are described by a reactor flux uncer- 
tainty <tr and two fiducial volume errors of the near and far o F detectors. Then Eq. (|Sj) 
will be replaced by 

# of bins 

X 2 (\a N ,a F ,a R ) = ^ ^ 2^(1 + a R + a d )T dti - O dji + 



i=l d=N,F 



+0 (iii log 7 — -i — — _ + ^ + io 

;1 + a R + a d )T dji J a R a z N ap 



1 For details on the event rate calculation, see Ref. [6]. 
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Figure 3: Luminosity scaling of the sin 2 2#i3 sensitivity for different systematical error assumptions (90% 
confidence level, figure similar to Ref. [19]). The thick curve corresponds to Eq. (TIT))) for ctr = <7 corr = 2.8%, 
ctjv = <Jf = funcorr = 0.6% and an additional calibration error. The figure is computed with example5.c 
provided with the GLoBES software. 

This function can be denned with user-defined systematics. Similarly, calibration errors, 
shape errors, uncorrelated bin-to-bin errors, etc., can be introduced. We show in Fig. [3] 
the result of such a calculation, where the thick curve corresponds to Eq. ffTOl) plus an 
additional energy calibration error|§ Furthermore, more complicated setups using multiple 
discrete sources and detectors can be simulated. For details, we refer to the manual [15]. 
In the past, user-defined systematics has (as experimental feature) been used in Ref. [19] to 
simulate Double Chooz, in Ref. [20] to simulate T2KK, and in Ref. [21] to simulate various 
reactor experiment setups and complicated backgrounds from geoneutrinos. 

In summary, GLoBES currently supports the following systematics \ 2 functions: 

chiSpectrumTilt Spectral information, signal and background normalization errors and 
tilts 

chiNoSysSpectrum No systematical errors, but spectral information 

chiTotalRatesTilt Total rates only, signal and background normalization errors and tilts 
chiSpectrumOnly Spectrum only (free signal normalization), no other systematics 

2 In fact, Fig. [3] was computed with Gaussian statistics for efficiency reasons, whereas GLoBES internally 
uses Poissonian statistics as a standard; see Ref. [15]. 
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chiNoSysTotalRates Total rates only, no systematics 

chiSpectrumCalib Spectral information, signal and background normalization errors and 
energy calibration errors 

chiZero Passive systematics (returns x 2 = Oj only provides the event rates for the access 
by other rules or experiments) 

Any other name User-defined systematics 

5 Adding external information with user-defined priors 

In order to include correlations and degeneracies, GLoBES provides (local) marginalization 
routines to project the fit manifold onto a subspace. For example, if one wants to compute 
the sin 2 2#i3 precision, one scans the sin 2 26^ direction and marginalizes for each fixed 
sin 2 26* 13 over the other oscillation parameters. This marginalization is performed after the 
systematics Xp U \\ m Eq. © (or similar) has been determined^] 

Xp roj ( A i) : = mm (Xp Ull (A 1; A 2 ) ) . (11) 

{A 2 } 

In this case, the parameter space of oscillation parameters A is projected onto the subspace 
Ai. If one wants to include external input, which only depends on the oscillation parameters 
and not on systematics, one can do so at this level. GLoBES provides the pre-defined 
possibility to include external Gaussian constraints on the oscillation parameters (and the 
matter density, which is treated as another oscillation parameter). For example, an external 
constraint on # 12 is added before the marginalization in Eq. (fTTj) by the replacement 

xuw - xuw + (gi2 ;/° 2)2 (12) 

a o 12 

with the central value 9\ 2 , where the input is added, and the la Gaussian error ag 12 . This 
corresponds to an external measurement of the solar mixing angle at the best-fit value 9\ 2 
with the error ag 12 . 

In practice, the external error on 6*12 may not be Gaussian or one may want to add other 
external information. Therefore, GLoBES 3.0 and higher provides the concept of user-defined 
priors /(A). Instead of Eq. (fT2l . we then have 

Xpuii(A) - XUW + /(A) • (13) 

One can easily imagine that this simple concept allows for a high degree of freedom. Ex- 
amples where this feature has been used as an experimental feature are the combination of 
terrestrial neutrino data with the information from neutrino telescopes [22] and the combina- 
tion of long-baseline data with atmospheric neutrino data [23,24]. Another very interesting 

3 The marginalization order systematics - correlations is, strictly speaking, not valid anymore for the new 
hybrid minimizer provided with GLoBES 3.0 as an experimental feature. However, we keep this description 
for pedagogical reasons. 
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Experiment 



File name 



Short description 



Refs. 



Superbeam experiments: 



T2K 
T2HK 

NOz/A 
SPL 



T2K.glb 
T2HK . gib 
NOvA . gib 
SPL. gib 



Reactor experiments: 



Reactor-I 

Reactor-II 

DoubleCHOOZ 

Beta beams: 
Low 7 
Medium 7 

Variable 7 



Reactorl . gib 
Reactor2 . gib 
D-Chooz_near . gib 
D-Chooz_f ar . gib 

BB_100.glb 
BB_350.glb 

BBvar.WC . gib 

BBvar.TASD . gib 



Neutrino factories: 



Standard 



Variable E u 



NFstandard . gib 
NFvar . gib 



Gold + Silver NF.GoldSilver . gib 
Hybrid detector NF_hR_lT . gib 



J-PARC to Super-Kamiokande [9,25] 

J-PARC to Hyper-Kamiokande [9, 25] 

Fermilab NuMI beamline off-axis [26, 27] 

CERNtoFrejus [24,28,29] 

Small reactor exp. , C = 400 1 GW yr [30] 

Large reactor exp., C = 8 000 1 GW yr [30] 

Double Chooz near detector [19] 
Double Chooz far detector 



7 = 100 CERN to Frejus scenario [24] 
7 = 350 "refurbished SPS" scenario (or [31] 
other accelerator) 

Variable 7 beta beam with water [32] 
Cherenkov detector, 50 < 7 < 500 
Variable 7 beta beam with totally active [32] 
scint. detector (TASD), 100 < 7 < 3 000 

Standard neutrino factory, 50 kt magne- [9] 

tized iron calorimeter, = 50 GeV 

Variable neutrino factory, disapp. chan- [9,33] 

nels without CID; 10 GeV < < 80 GeV 

As NFvar . gib plus 5 kt ECC detector for [9, 33, 34] 

Silver Channel measurement 

As NFvar. gib, but lower threshold and [9,33] 
better energy resolution 



Table 1: Pre-defined experiment prototypes coming with GLoBES 3.0: The columns represent the AEDL 
filename, a short description, and the references in which the files are originally used and discussed (except 
for minor modifications, such as a different implementation of the energy threshold function). 

option is the use of penalties in degeneracy localization, as suggested by Thomas Schwetz. 
For example, if the minimizer runs into the wrong octant, user-defined priors can be used 
to add a penalty and prevent it from doing so. 

6 AEDL changes and experiment prototypes 

AEDL (Abstract Experiment Definition Language) is a special syntax developed for the 
description of experiments in GLoBES, represented by simple text files. Core part of AEDL 
are the following constructions (for details, see Ref. [15]): 

Channel A channel represents an oscillation channel from the source flux, over the initial 
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and final flavors and the use of neutrinos or antineutrinos, to the interaction type and 
cross section. The result is an event rate vector. 

Rule A rule combines one or more signal channels with one or more background channels. 
The rates from these channels are added, and the signal and background are imple- 
mented with specific systematics, i.e., systematics is connected to the concept of a 
rule. The result of a rule is a systematics x 2 corresponding to Eq. ((7j). 

Experiment An experiment is a combination of one or more rules, which may correspond 
to different appearance and disappearance channels, neutrino and antineutrino oper- 
ations modes, etc.. An experiment has a common detector and baseline, but may 
use several different fluxes (such as for neutrinos and antineutrinos). Therefore, the 
matter density profile is the same for all rules of an experiment. The result of an 
experiment is again a systematics x 2 ■, but summed over all rules. 

Experiment combination An experiment combination consists of different experiments 
with different detectors, baselines, and (uncorrelated) matter density profiles. The 
result of an experiment combination is again a systematics x 2 , but summed over all 
experiments. 

Note that after the x 2 is summed over all rules or experiments as specified by the user in the 
application software, Eq. (Hip is applied, and the oscillation parameters are marginalized 
over. 

In GLoBES 3.0, this construction of experiments is extended by user-defined systematics, 
which allows the user to break up the strict coupling between systematics and rules. For 
example, systematics can be correlated between different rules or experiments, such as for 
a reactor experiment with identical near and far detectors. In addition, one can introduce 
user-defined systematical errors. For that purpose, the former error dimension concept has 
been changed and extended. For example, one has to give each user-defined systematics a 
name, and this name has to be matched by the application software. Further changes in 
AEDL are the requirement to specify the minimum GLoBES version the AEDL file can be 
used with, and the extensions by lists as AEDL variables and by an interpolation feature. 

As far as the experiment prototypes are concerned, a list of the ones provided with GLoBES 3.0 
is given in Table [TJ These prototypes can be modified for the definition of one's own 
experiments (experimentalist), or they can be used for the test of new physics concepts 
(phenomenologist /theorist). Note that some of these experiments may contain integrated 
luminosities, baseline, fluxes, efficiencies, or other factors, which the user may not agree with 
and which can be easily adjusted. Some of the AEDL files have been changed or updated 
from earlier versions, such as for T2K and T2HK, or the standard neutrino factory. Other 
files represent the results from recent developments, such as the variable energy beta beam 
or neutrino factory files. They push AEDL to the edge in terms of advanced features, such 
as variable bin widths, use of the new interpolation feature, threshold function implemen- 
tations, etc.. In particular, note that Double Chooz requires user-defined systematics, and 
the beta beams require built-in beta beam fluxes. Both of these features are only available 
in GLoBES 3.0 and higher. Some of the files require that AEDL variables be pre-set by the 
user, such as the muon energy for the variable energy neutrino factory. These variables are 
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described in the corresponding file headers. Old AEDL prototypes should still run, but will 
not be maintained by the GLoBES team anymore. 

7 Internal changes and improved performance 

Sophisticated GLoBES simulations with several experiments and complicated parameter 
correlations can take several days to finish even on modern computers. To mitigate this 
problem as far as possible, we have taken several steps to optimize the performance of the 
multi-dimensional minimizations, which are at the heart of each GLoBES x 2 analysis. 

Firstly, the GLoBES minimizers distinguish between minimization over oscillation parame- 
ters and minimization over systematics (nuisance) parameters. This is necessary because 
each step in the oscillation parameter space requires a time-consuming re-computation of 
oscillation probabilities, while a modification of the systematics parameters does not. The 
standard procedure is to perform a full minimization in systematics space after each step in 
oscillation space. This avoids interference of the two categories of parameters, and is known 
to have excellent convergence properties. The numerical method that does the minimization 
is the Powell algorithm [35] . 

In GLoBES 3.0, we provide an alternative method, which does the full minimization in one 
go, and is therefore several times faster in most cases. In this algorithm, each Powell iteration 
over oscillation parameters is followed by one iteration over systematics parameters. Since 
there may be some pathological cases where this interleaved minimization strategy behaves 
differently than the old method, it is by default not used in GLoBES 3.0. This ensures 
maximal compatibility with old application programs, but it requires users who wish to 
benefit from the new algorithm to explicitly select it by an API function call. 

Besides optimizing the minimization functions, it is also mandatory to optimize the routines 
which are called by them to calculate x 2 f° r one specific set of parameters. In particular, 
the computation of oscillation probabilities, which has to be done every time the oscillation 
parameters have changed, should be as efficient as possible. The most expensive step here 
is the diagonalization of the hermitian 3x3 neutrino Hamilton operator in matter. The 
LAPACK routines [36], which were used to solve this eigenproblem in previous versions 
of GLoBES, turned out not to be the optimal choice because they are mainly optimized 
for very large matrices [37]. In the case of the small and usually well-conditioned matrices 
appearing in GLoBES, they spend too much time evaluating their parameter list and planning 
the optimal diagonalization strategy. Therefore, GLoBES now relies an algorithm which has 
been designed specifically for the diagonalization of hermitian 3x3 matrices. It is essentially 
the extremely fast analytical algorithm discussed in Ref. [37], but for the rare cases where 
this may become inaccurate, it falls back to the QL method [38,39]. 

As far as the installation of GLoBES on different systems is concerned, GLoBES has been 
improved in the past, such as to handle different Linux distributions. In addition, for high- 
throughput computing, GLoBES has been used or tested on 64 Bit systems, Condor clusters, 
and parallel clusters. An important prerequisite for the parallelization is the option to build 
static binaries, which is now illustrated in the examples. With these improvements, the 
performance of GLoBES has been pushed to the edge. Ref. [33] is probably the currently 
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most demanding product of GLoBES, for which most of these features have been developed. 
8 Summary and conclusions 

We have presented GLoBES 3.0, which is a simulation software for future terrestrial neu- 
trino oscillation facilities. GLoBES is now no further restricted to one source-one detector 
configurations, but it can simulate the systematics of arbitrary combinations of discrete sets 
of sources and detectors. In addition, user-defined systematical errors can be introduced. 
The new concept of user-defined priors allows the inclusion of external information at a 
different level. For example, constraints from atmospheric neutrino experiments or neu- 
trino telescopes can be taken into account. In addition, they provide a powerful handle on 
complicated degeneracy localization. On the probability level, GLoBES now supports the 
simulation of non-standard physics beyond three-flavor neutrino oscillations without a re- 
compilation of the software. A very important prerequisite for this feature is the ability of 
GLoBES 3.0 to carry more than six oscillation parameters. New experiment prototypes make 
extensive use of these new features, as well as they are illustrated by examples in the man- 
ual. One prominent example is the Double Chooz experiment, which requires user-defined 
systematics to simulate correlated errors between the near and far detectors. 

In conclusion, GLoBES 3.0 provides a number of new improvements, such as by including 
more flexibility. At the end, the developers of a software can never exactly foresee what the 
users are up for, and therefore flexible concepts allow for more space in problem realizations. 
In this version of GLoBES, three new such flexible concepts have been introduced: user- 
defined oscillation probabilities, user-defined systematics, and user-defined priors. These 
approaches will make an easy adjustment of the software possible for both experimentalists 
and theorists, and will allow for flexibility at all levels. Finally, note that the past has shown 
that there is often more than one possibility how a specific problem can be implemented in 
GLoBES. 
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